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ABSTRACT 


The theory of inelastic scattering is formulated, using detailed nuclear wave func- 
tions which are obtained from nuclear structure calculations. The inelastic scattering 

is assumed to proceed through a "direct” reaction, and nuclear states are obtained by 
o 

projecting states of good (J , J z ) from a variational Hartree-Fock wave function. Gen- 
eral expressions for transition amplitudes and the nuclear form factors are derived, and 
the form factors for inelastic proton scattering are presented. These form factors are 
compared with those based on the shell model and macroscopic collective model. It is 
concluded that the restrictions imposed in the Hartree-Fock model must be relaxed if the 
model is to give the correct order of magnitude for reaction cross sections. 
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INELASTIC SCATTERING FORM FACTORS USING PROJECTED 
HARTREE-FOCK WAVE FUNCTIONS 
by Richard C. Braley and William F. Ford 
Lewis Research Center 

SUMMARY 

The theory of inelastic scattering is formulated by use of detailed nuclear wave func- 
tions which are obtained from nuclear structure calculations. The inelastic scattering is 
assumed to proceed through a "direct" reaction, and nuclear states are obtained by pro- 

9 

jecting states of good (J , J ) from a variational Hartree-Fock wave function. 

z 

General expressions for transition amplitudes and the nuclear form factors are de- 
rived, and the form factors for inelastic proton scattering are presented. These form 
factors are compared with those based on the shell model and macroscopic collective 
model. It is concluded that the restrictions imposed in the Hartree-Fock model must be 
relaxed if the model used here is to give the correct order of magnitude for reaction 
cross sections. 


INTRODUCTION 

The success of the distorted wave method in direct reactions, coupled with the avail- 
ability of more accurate data, makes possible detailed studies of the effects of nuclear 
structure in nuclear reactions. In the past, however, only relatively simple nuclear 
wave functions have been used in such studies, namely, those based on the spherical- 
shell model or the macroscopic collective model. The use of these models in reaction 
studies has been discussed exhaustively in the literature (refs. 1 to 3). Although the re- 
sults so obtained have been fairly successful, it is well known that such primitive nuclear 
wave functions are inadequate for nuclear structure calculations. This raises consider- 
able doubt as to their reliability for the scattering calculations. The point has been em- 
phasized by Glendenning and Veneroni, who, in a recent study of inelastic proton scatter- 
ing from the even nickel isotopes, found that inelastic scattering is rather sensitive to 
the details of a microscopic description of the target (ref. 4). Satchler (ref. 5) has given 
a discussion of inelastic scattering based on the shell-model description of nuclei situated 



near closed shells. Reaction studies (ref. 6) which were done after Satchler’s work, 
namely, A(pp')A* and A(p, n)B, indicate that a more detailed description of the nuclear 
form factor is required. 

In reaction theory, efforts to use detailed descriptions of states of a collective nature 
have been minimal (ref. 7). Furthermore, it is well known that the spherical shell model 
does not produce collective effects. Detailed descriptions of the inelastic scattering of 
a -particles and deuterons through collective excitation are generally limited to vibra- 
tional states in the region of spherical nuclei. At the present stage of development, re- 
action studies which make use of the collective model provide only the deformation pa- 
rameter and multipolarity of the transition. In nucleon inelastic scattering, if such a 
model is used, there is no provision for the study of spin exchange or particle exchange, 
and the absence of spin-dependent forces prohibits the study of unnatural parity states 
which may be excited via spin-flip. More complex projectiles, such as a-particles, may 
require the use of coupled-channel methods to study parity states. 

When the macroscopic model is used in distorted wave Born approximation (DWBA) 
reaction calculations, it is assumed that the inelastic processes result from nonspherical 
components in the optical potential. A Taylor expansion of the optical potential, in terms 
of the deformation parameter, is made, and the first order term in the expansion is re- 
tained as the interaction which gives rise to inelastic scattering. In the microscopic de- 
scription, the choice of effective interaction is not this simple. A sum of effective two- 
body forces between the projectile and the target nucleon will be assumed in this paper. 
The justification of such a choice will only be apparent after the results of detailed calcu- 
lations are compared with experiment. As Satchler has suggested (ref. 5), the use of 
such an effective interaction assumes that multiple scattering is not important; however, 
if one uses phenomenological optical potentials to describe the distorted waves, then 
multiple scattering is neglected only in the off-diagonal parts of the scattering matrix. 
There is little that can be said about the analytic form of the effective interaction because 
it cannot, as yet, be calculated. The nucleon-nucleon interaction may be chosen to be 
either the two-nucleon t-matrix or a spin -dependent phenomenological interaction. The 
nucleon-alpha or nucleon -deuteron t-matrix may be chosen in a similar fashion. It is 
suggested, however, that for the case of nucleon-nucleon scattering the phenomenological 
projectile -target nucleon interaction (for intermediate energies, 10 to 50 MeV) should not 
be radically different from the residual two-body force in the structure problem because 
the interaction occurs in the field of nuclear matter. 

In recent years, the Hartree-Fock (HF) theory (ref. 8) has received considerable 
attention in nuclear structure studies in the (p, s-d) shell nuclei. Recently Bouten, et al., 
using HF methods, were quite successful in a structure study of carbon 12 (ref. 9). 
Structure calculations for other nuclei (refs. 10 and 11) in the (1, s-d) shell have been 
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equally successful. The self-consistent HF method provides a microscopic model which 
may be used to study systems of many bound nucleons. A particularly attractive feature 
of this model is the inclusion of residual two-body forces between the nucleons in the 
system. The self-consistency requirement coupled with the use of residual two-body 
forces in the Hamiltonian for the nucleus gives rise to deformed orbits, if the angular 
HF method is used. If a static (time independent) HF calculation is performed, one ob- 
tains an intrinsic wave function for a statically deformed nucleus, and it is possible to 
project from the intrinsic wave function states of good J and M which exhibit a rota- 
tional character. The extension of these methods to a time dependent formalism may be 
used to describe vibrational states (ref. 12), although these states are not discussed in 
any detail herein. 

The major objective of this paper is to formulate a microscopic theory of inelastic 
scattering from deformed nuclei. The nuclear wave functions for the target and residual 
nucleus are obtained by projecting states of good J and M from the intrinsic HF state. 
These wave functions are then used to calculate nuclear form factors for reaction studies. 
Furthermore, a recent study of the exchange term in the transition amplitude for inelas- 
tic nucleon scattering indicates that the exchange term is not negligible (ref. 3); these 
effects will be included in the reaction formalism to be presented herein. 

The transition amplitude is discussed in detail, and nuclear form factors are defined; 
detailed expressions are given for the nuclear form factors for the case in which the nu- 
clear states are described by projected HF wave functions. A comparison of the HF 
form factors with those obtained from use of projected shell model states and the macro- 
scopic collective model is also presented. 


SYMBOLS 


A(B) refers to (A-l), particle systems 

A. ,B. irreducible tensors of rank j 
jm’ jm J 


a,b, . . . 

a s 

B 

B XjU 


C A 

jm 


quantum numbers for particular single -particle states 
diffuseness of Wood-Saxon well 
N x N matrix 
element of matrix B 

normalization factor for macroscopic form factor 
expansion coefficients for Hartree-Fock orbits 
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“mkW 

e A 

f ab (r ) 
LS V o' 

h 

i(f) 

(J,M) 

(JpM.) 

(ja m a> 


rotation matrix 

Hartree-Fock single -particle energy 

single -particle form factors 

Hartree-Fock Hamiltonian 
incident (outgoing) channel 

total angular momentum transfer and its projection onto z-axis 
total angular momentum and projection for nucleus in channel i 

total angular momentum and projection for single-particle 

Clebsch-Gordan coefficient 


K 

k i 

(lm l ) 

Vlh 

«L S j| 

l 2 S 2h 

v- J 
m 

/ 

P L (r 0' F P 





projection of J onto z-axis of the body-fixed system 
wave number of the projectile in channel i 
orbital angular momentum transfer and projection 
single -particle orbital angular momentum and projection 

symmetric 9-j coefficient 

reduced mass of the projectile 
normalization factor 

Legendre polynomial with argument which is cosine of angle between 
r Q and r ^ 

angular momentum projection operator 
projection operator for singlet-odd states 

projection operator for triplet-even states 

Heisenberg exchange operator 

radius of Wood-Saxon well 

spectroscopic factor for angular momentum transfer J 
neutron part of spectroscopic factor 
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(s, Mg) 

(SqM) 


t jm 

qr 

LSJM t 

V 

V 0 * 

t™(0) 
t(0, 1) 

u 0 l 3 2 j 34; 33 J') 


proton part of spectroscopic factor 

spin angular momentum and projection 
projectile spin and projection onto z-axis 

single -particle transition amplitude 
direct -reaction transition amplitude 

coupled form of transition amplitude 

transition amplitude for transfer of definite L and S 
effective single -particle interaction 

coupled form of t^ r (0) 

projectile-nucleon interaction 

[(2j 3 + 1)(2J’ + 1)] 1/2 J j 4 ;j 3 J’) 


Uq strength of real part of optical potential 

V 0 p(r^) optical potential 

v(i, j) two-particle interaction used in structure calculations 

W(jfj 2 J j 4 ;j 3 j') Racah coefficient 


Y LM L ^ r 0^ 




M 

LSJ 




X 


CT S 

dcr/d£2 

T 

<3> or 


spherical harmonic with argument denoted by unit vector Tq 

spin-angle tensor 
deformation parameter 

destruction and creation operators for single -particle states 

Hartree-Fock single -particle orbits 
spin tensor operator 
differential cross section 
iso spin 

determinantal wave function for A-particle system 



r-J 


$ or 

*K 



*(1,. 

. ,,A) 

*ia< 2 ’ 

• • • > A) 




Si 


rotated Slater determinant 
single-particle basis function for particle 1 
Hartree-Fock orbit 

rotated Hartree-Fock orbit 

nuclear wave function for nucleus containing A nucleons 
nuclear wave function in channel i with particle 1 

distorted waves for incident (+) and outgoing (-) channels 

Euler angles 

first two-plus state 


TRANSITION AMPLITUDE 
General Formulation 

Under the assumption that the projectile-target interaction can be written as a sum 
of two -body interactions, the DWBA transition amplitude for inelastic scattering from a 
target containing A nucleons is 

Tjf ~ A^f"^ (0)^(1 • • • A)|t(0,l)|^ +) (0)^.(1 . . . A)) (1) 

The initial and final nuclear states are represented by T. and 'f^, and the initial and 

final projectile states by the distorted waves and We assume that integra- 

tion over the internal coordinates of the projectile has already been carried out, so that 
0 stands for the spin, isospin, and center-of-mass coordinates of the projectile. The 
projectile-nucleon interaction is represented by t(0, 1); it may be the t-matrix, for in- 
stance, or a phenomenological pseudopotential. For nucleon-nucleus scattering, the 
major exchange term (ref. 3) may be included simply by writing 

t (0, l)(l - Pqj) 

XJ 

in place of t(0, 1), where Pqj is the (Heisenberg) operator which exchanges all the co- 
ordinates of particles 0 and 1. 
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Because of the assumption that the projectile -target interaction is two body in char- 
acter, it is possible to express the transition amplitude in terms of single -particle tran- 
sition amplitudes. We begin by expanding the initial and final nuclear states in a complete 
set {<Pj m } of single -particle wave functions, where j stands for all the quantum num- 
bers (nljT) except for the z-component of angular momentum, which is denoted by mg. 
(Often, however, a will be written for the set j m , the set j , or even the set n„l t„, 
when no confusion will result. ) Thus, 

*i<l • • • A) = A' 1 / 2 Yj <P a (1 )*ia (2 • • • A ) (2a) 

a 

* f (l . . . A) = A'' 1/2 Y, ^b (1 )% (2 • • • A ) ( 2b ) 

b 

The scattering amplitude then becomes 

T if = S<*£blV T ab < 3 > 

ab 

where the single -particle transition amplitude T ab is given by 

T ab = <Vf -) (0)^ b (l) |t(0, l)|^W a (l)) (4) 

It is convenient to express both Ty, and T ab as sums of amplitudes for transfer of total 
angular momentum (J,M)'*' 


JM|j f M f )Tj M 

JM 

(5a) 

£<ia m a.™lvV T ™ 

JM 

(5b) 


^Condon and Shortley conventions are used for the Clebsch-Gordon coefficients 
(Jp nip ] 2 m 2 I JM> . Reduced matrix elements, W-coefficients and the 9-J coefficients, 
are defined in ref. 13. The U-coefficient was introduced by H. A. Jahn (ref. 14). 
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And also to express \f^ a and in terms of angular momentum eigenstates: 




, j a m a 


j A m A 


(6a) 




fb 


j b m b 


(6b) 


Then the relation 


= Z s J (lf i ab > T ™ 
ab 


(7) 


holds, where the spectroscopic amplitude Sj is given by 



( 8 ) 


In the event that the distorted waves do not have spin-orbit coupling, one can decompose 
still further into amplitudes for transfer of definite L and S: 


tm , , LSJM t 

T diVi = > (S Q M f , SMg] S Q M i ) (LM l , SM g j JM)T u (9) 

lm l 

SM s 


LSJM t lsjm t 

t u -2 J Sj<u|ab,Tab <10) 

ab 

Here Sq is the spin of the projectile, and and are its initial and final projec- 
tions on the z-axis. The sum over initial and final states of the entire system is inco- 
herent in this case, so that we obtain 
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( 11 ) 


where 


'da\ = / m \ 2 h 2J f + 1 g LSJ 

dfi /if " \27rR 2 / k i 2J i + 1 (2L + 1K2S + 1} 


a LSJ 


Z LSJM 
T “ 


2 


M 


L 


(12) 


Single-Particle Transition Amplitude 

With equation (7) the first step of the program, namely, to separate details of the 
reaction mechanism (given by T ab ) from details of the nuclear structure (given by Sj), 
is complete. The single -particle transition amplitude is of secondary concern; it may 
be evaluated in a number of ways (ref. 6), depending on the specific forms of t(0, 1) and 
the distorted waves. For example, when the optical potential has no spin-orbit term, 
and the spin dependence of the two-body interaction can be expressed in terms of the spin 
tensors Oq = 1 and = 2S as 


t(0, 1) t g (r 0 , r^<7 s (0) • CTg(l) 
S 


the equation 


LSJM, 

r ab 


“ ( S 0 a S I S o) g LS,I T ab 


lsm t 


is obtained, where 


gLSj = ^2(2^ + l)(2j a + 1)(2S + 1)(2J + 1) 



!/ 2 0 

S J f 

V 2 lb. 


(13) 


(14) 


(15) 


and 
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LSM 
^ab 


L 2L 


2 *b 


2E 


;'a m i a ' LM L 


l h m U 


I ‘s(V F l) | 


m 


V”*b 


( 16 ) 


If tg is central, an expansion in Legendre polynomials 


t s( r 0’ r l) + 1 P L( r 0 ‘ r l) t LS^ r 0’ r l^ 


yields 


LSM t 

r ab 


- /4 \? 0 ) f LS^ r O^ Y LM L ( r o)^i + \ r o) dr O 


ab 


where the single -particle form factor f Lg is given by 


(17) 


(18) 


f LS^ r 0^ _ Y lHV f ^b^ r l^LS^ r 0’ r l^a^ r l^ r l dr l ( 19 ) 


A convenient way to express this result is to define a nuclear form factor, 
such that 


f if ( r ) 
LSJ' 0'’ 


T if L “ ^ S oU CT sll S 0^/^f ^( r o) f LSJ^ r O^ Y LM L ( r o) 1 ^ / i + ^( r o) dr 0 ( 20 ) 

and then to write it in terms of the single -particle form factors: 

f LSj( r o) = S S J^l ab ) g LSJ f Ls( r o) ( 21 ) 

ab 

If the two-body interaction t(0, 1) is isospin-dependent, t^g^, r^J will be a function 
of the isospin coordinates of the projectile and the bound nucleon. One could introduce a 
general formalism for handling the isospin at this point. When Hartree-Fock wave func- 
tions are used, however, there occurs a natural separation of the spectroscopic amplitude 
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Sj into a neutron part Sj and a proton part S^. Examples of the isospin treatment will 
therefore be deferred to the section Spectroscopic Factor. 

Exchange term for inelastic nucleon scattering . - In the case of inelastic nucleon 
scattering, Amos, et al. (ref. 3) have shown that the exchange term is nearly in phase 
with the direct term, and makes a substantial contribution to the cross section. To in- 
clude this term, when the spin-dependence of t(0, 1) is as in equation (13), tg^Q, r^ in 
equation (16) must be replaced by 



MT 

where P^ is an operator which exchanges the space and isospin coordinates of nu- 
cleons 0 and 1 and where 


e S o. = (-1) S+S, " 1 (2S’ + l)w/— III 

\2 2 2 2 


SS' 


_ 1 2 . . 

= si s’l 


(23) 


The exchange terms cannot be put in a form like equation (18), of course, so no form 
factor can be defined. 

Form factors for coupled channel calculation . - In some reactions, transition rates 
to various excited states of the nucleus are comparable in magnitude to the ground state 
(elastic) transition rate. Consequently, even an approximate wave function should have 
considerable projection onto these excited states, and the resulting analysis is a coupled- 
channels calculation. For such a calculation a set of effective projectile -target interac- 
tions of the form 


t^O) = A ('& f (1 . . . A)jt(0, ... A)) 

= Yj < J i M p JM I J f M f >4 M (0) (24) 

JM 

is needed, where the integration extends over the nuclear coordinates and the internal 
coordinates of the projectile and i and f now label the nuclear states included in the 
expansion of the total wave function. Proceeding exactly as before, a single -particle ef- 
fective interaction is introduced. 
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tab* 0 ) s ^b (1) l t(0 ’ 


= Z < j a m a» JM |jb m b )t 2 I(0 ) (25) 

JM 


The two equations are related by 

4 M (°) =Z Sj^ (26) 

ab 

When t(0, 1) can be expanded in spin tensors and spherical harmonics as in equations (13) 
and (17), the equation 


4 M (°) = Z (- 1 ) L+S+Jf ^SJ (r o)(- 1 ) M ^LSj(°) (27) 

LS 

is obtained, where f^SJ is de * ined in equation (21) and the spin-angle tensor is 

given by 

Yj < LM L' SM sl JM > Y LM L < r 0> < 'SM s (°) < 28 > 


Spectroscopic Factor 

The evaluation of Sj is performed using Hartree-Fock wave functions. As is well 
known, the Hartree-Fock method consists of approximating the nuclear wave function by 
a single Slater determinant, 


$= (A!)" 1 / 2 det{<p A } (29) 

where {cp^ } is a set of single -particle wave functions, or orbitals. These functions 
are determined by the requirement that <3> minimizes the expectation value of the nu- 
clear Hamiltonian H, which leads to the coupled eigenvalue equations 


h( P\ = e \<P\ 


(30) 
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where 


h(i) = ^ + X Art®* “(* ' p «) , v (i) (3i> 

M 

with the M-sum restricted to those orbitals used to construct <$. 

The intrinsic HF wave function does not depend on the orientation of the deformed 
field relative to the space-fixed system, and, clearly, the single -particle energies of the 
nucleons in the deformed field do not depend on this orientation. Thus, states of the in- 
trinsic system which are based on the same deformed field are degenerate with respect 
to orientation; any two such states are simply related by 

<P K (x’) = R( fi )<P K (x") (32) 


where is the orientation of x M relative to x'. Also, it should be noted that the 
states <p K are not states of good J. However, it is possible to construct states of good 
J and M in the space -fixed system by taking appropriate linear combinations of the de- 
generate states <p K (x'). 

In order to determine the expansion coefficients, a particular <p K (x') is expanded in 
a complete set of states ji//^(x')j which are normalized eigenstates of and J z in 
the body-fixed system: 


<P 


K 




( 33 ) 


K, 


Here. x f stands for all the coordinates in the body -fixed frame. If tf/v* (x T ) is now ex 

TV/r J 


pressed in terms of the states i^(x) in the space-fixed frame 


^j(x’) = R(Si)*J(x) = 2 D ^ K (Wf (x) 

M 


( 34 ) 


where £2 is the orientation of the body-fixed frame relative to the space-fixed frame, 
then 


JM 


( 35 ) 
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This expression may be inverted to give 


/»< k ( 8 )V) (36) 

8 ” 2 “j 

A linear combination of intrinsic states, with 

D MK (B) * (37) 

8 " 2 “j 

as the weighting coefficients, gives rise to a state of good J and M in the space-fixed 
frame (eq. (17)). It is not necessary to calculate each cp-^ix') individually, because 

<P K (x') = R (*% K (x) 


hence, equation (17) can be written 

^’(3) = f JSiD^ K (Si)’R(fi)<,, K (3) (38) 

81,2 “j 

Thus, it is sufficient to calculate cp-^ in the space-fixed frame and carry out integration 
(18), which is denoted symbolically by 

^(x) = PmK^K® (39) 

It is clear that, from one intrinsic state <p^, many states may be generated; these 
are usually associated with the rotational bands found in deformed nuclei. The Slater de- 
terminant <p£-(x) may be written as the antisymmetrized product of two determinants, 
one of which represents the spherical closed-shell core and the other the M nucleons 
outside the core. Because the core determinant has J = 0, it will be unaffected by the 
rotation operator R(£2). Hence, the net effect of the rotation operator in equation (18) 
will be to rotate each of the orbits in the extracore determinant only, or equivalently, 
to rotate only the extracore orbitals in cp £. 

Of course, the state 'T'j M yields a higher expectation value of H than $ does, be- 
cause of the cross terms, but it is hoped that the discrepancy is not large. Calculations 
in the s-d shell confirm this hope, at least for low-lying nuclear levels (refs. 10 and 11). 
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Abetter procedure would be to determine the set { <p^ } by minimizing <^|h|^) di- 
rectly, but this is a formidable task and until recently (ref. 9) has not been attempted. 

Quite commonly the nucleus is represented as N nucleons outside a closed-shell 
core, and the core orbitals are not varied in the Hartree-Fock procedure. In that case 
equation (37) becomes 

h(i) = h c (i) + ^ f dT.<p+(j)v(i, j)^l - P^<P M (j) (40) 

M 

where h c (i) is the Hamiltonian (eq (31)) for the core orbitals and the sum over p is 

limited to the extracore orbitals in 4>. Usually h £ is not dealt with directly; instead, 

its eigenfunctions are assumed to be shell model orbitals (harmonic oscillator functions, 

for instance), and its eigenvalues are taken from experiment (ref. 15). As before, eigen- 
o 

states of J and J z are obtained using equation (32), but, because of the closed- shell 
core, a simplification occurs. The determinant $ may be written as the antisymme- 
trized product of two determinants, one of which represents the core nucleons. This 
latter determinant has J =0 and is thus invariant under rotations. Hence, the net ef- 
fect of the rotation operator in equation (32) is to rotate only the extracore orbitals in <3>. 

Whatever method is used to determine {cp^ } , we assume that the nuclear wave func- 
tions are single Slater determinants projected onto states of good J and M by means of 
equation (18). We also assume the existence of a closed-shell core, identical in initial 
and final states, because, in the final results, the core can always be taken to contain no 
particles. Finally, we assume that the set {<p^} has been made orthonormal, because 
any nonorthogonal components would make no contribution to the Slater determinants. 

The first task is to show that the summation over single -particles in equation (3), 
and consequently in equation (7), is restricted to states outside the core. This is easiest 
to show in the notation of second quantization, using fermion creation and annihilation op- 
erators, for, then, the overlap integral in equation (3) may be written 

<%bKa> = <%l^al'V 

= 6 ab <*f I *i> - <*f I V*b I V (41) 

The first term on the right side of equation (34) vanishes for inelastic transitions, be- 
cause the initial and final nuclear states are orthogonal. The second term will also 
vanish, if either a or b refers to a core state, because the core has been unaffected 
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by the projection procedures of equation (32) (except for antisymmetrization) and is con- 
sequently still filled. 

We next turn to evaluation of the spectroscopic amplitude Sj. If the nucleus is not 
spherical, the orbits cp ^ will, in general, be of the form 

n = < 42 > 

jm 

where { cp^ m } is the complete set used earlier. Consequently, an expansion of the de- 
terminant <$ leads to 


* = A- 1/2 £ C 4 m »> jm (l)* x (-1) 


Ajm 


(43) 


where 4>^(-l) is the cofactor of cp ^ in $ and is itself a normalized Slater determinant 
for A - 1 nucleons. 

The projection operator must now be applied to $. From the definition of 4>, 

it may be shown that P^^ picks out states whose and J z eigenvalues are J(JH) 

and K, respectively, and then changes the J eigenvalue from K to M. Consequently, 

J z 

Pmk transforms under rotations like an irreducible tensor Aj-^. The product of an- 
other irreducible tensor A. with PjLt^t may therefore be written 

jm ivi jv 


A jm P M'K’ = Yj 


J 1 M 1 


and application of P^-^ then yields 


(44) 


J’ 


P k A im P M-K- = (J'M 1 , jm [ JK) Bjm ( j; J'K t ) 


(j’M*, jm j JK) {JVjm’IJMjA^.pJ^, 
m 'p 


(45) 


(Note that Bj^ has been singled out and then transformed to Bjj^. ) Setting K r = M* 
and summing both sides of equation (45) over all J' and M' results in 
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(46) 


P MK A jm=£ A jm- £ (J'M'.jm | JK> <JV, ]m’| JM> pJ' M , 


m' J'M’m 


which is the desired theorem. 

Using equations (43) and (46) to construct the nuclear wave function yields 

^JM s 2 P MK $ 

K 

= A' 1/2 £<?.(!) Yj (JV,im|jM) Yj (47) 


3 m 

jm J|i 


u jm 1 
XKm’M’ 


From the definitions of 4> ia and , given by equations (2) and (6), it follows 


by comparison with equation (47) that 


/a Ma 

ia 




x l M A m a K i 


% 


j b m b 


fb 


Yj C b 2 < J B M B>jb I "bl J f K t >P M B M^*fX 2 (- 1 ) < 48b > 


A 2 M B m b K f 


, 'A t °A\ 
V P fb *ia ) 


The overlap integral needed for evaluation of Sj is therefore given by 

= X ( C b 2 ) C a 1 X <J A M A^a m al J i K i> 

M A M B 

K iK f 

X ibmblJfKf) (%x 


X, X„ m„m, 
12 a b 


M B M A 




(49) 


For simplicity, the normalization of the initial and final nuclear states has been ig- 
nored in the derivation of equation (49). It is clear, however, that 
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will not be normalized to unity, even if is, and equation (49) should therefore be mul- 
tiplied by the normalization factor 


where 

>i= = 

Jff = <*f |*f> = 


If 


KK’ 


I 

KK’ 




KK' 


KK' 






(50) 


(51a) 


(51b) 


The integrals appearing in equations (49) and (51) are evaluated by means of the inte- 
gral representation (32) for Consider the normalization integral, which contains 

terms of the form 

f d«Dj IK (n)*<$|i> (52) 

where 3 = R(fi)4>. The effect of R(£2) on <3> is to rotate each of the N extracore orbi- 
tals in 4>. The core orbitals in 4> and $ are identical, and so we find that 

<$|$)=det{B} (53) 

where B is an N x N matrix with elements 

B AM = <^1 V = < 54 ) 

The integral in equation (49) is treated in a similar fashion, except that now terms 
of the form 




% 




( 55 ) 
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are encountered. Here and <1>^ are normalized Slater determinants for A - 1 par- 
ticles, again with identical core orbitals. It can be shown in this case that 
is equal to the cofactor of B^ in the determinant det{B}. An elementary theorem 
concerning cofactors of determinants then permits us to write 


= V det{B} 
-1 


(56) 


where b ^ is an element of the matrix b = B 

Most Hartree-Fock calculations are made without isospin mixing, which means that 

each orbital is either a pure neutron state or proton state. The summation in equation 

(49) is therefore restricted by the condition t = t , and the result vanishes unless 

X 1 x 2 


a t^. One may thus speak of the proton part Sj or the neutron part Sj of the spec - 


T_ = 


troscopic factor defined in equation (8). 


FORM FACTORS FOR 2s- Id SHELL NUCLEI 
Structure of s-d Shell Nuclei 

o 

The general procedure that is used to obtain nuclear states of good J and was 

discussed in the previous section. In the calculation of the "intrinsic" state defined in 

equation (29), we have assumed that oxygen 16 constitutes the inert closed-shell core. 

The 2s- Id shell nuclei are considered to be axially symmetric, although there have been 

24 S2 

some calculations which indicate regions of asymmetry (namely, Mg , S ). However, 
because the prediction of asymmetric regions in the s-d shell is based only on the exami- 
nation of the intrinsic spectra and because the assumption of axial symmetry has already 
yielded excellent structure results (ref. 8), we believe that the assumption is justified. 

In this case, equation (32) becomes 


^JM - P MK $ K 

and the orbitals will be of the form 

i 


(57) 


(58) 


The subspace of the basis states is limited to the 2s-ld shell (i. e. , ldg^ ^ s i/2’ 
ldg/g), with harmonic oscillator radial dependence. The spherical single -particle 
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energies are taken from experiment, and the two-body interaction which appears in equa- 

o 

tion (33) is chosen to be a Rosenfeld mixture with a Gaussian radial dependence 

v(i, j) = V 0 e"^ r2 i (0. 3 + 0. 7 Oj - a.) 

3 

with Vq = +55 MeV. The harmonic oscillator constant, u = mco/R, was chosen such that 
\i = v. There has been no correction for center-of-mass motion. 

The determinate wave function <f> K is constructed in such a way that, along with 
each single -particle state, one also includes the time reversed state 

<59) 

j 

that is, a proton and neutron are put in the state (58) and the state (59). This gives rise 
to a fourfold degeneracy of each orbit in an even-even nucleus; the resulting determinate 
wave function is then the time reversed invariant. The construction of such states for 
odd nuclei is more involved and will not be considered. 


TABLE I. - ENERGIES AND WAVE FUNCTIONS FOR NEON 20 
AND MAGNESIUM 24 


Nucleus 

Hartree-Fock 

Projection 

Projection of Hartree- 

Fock orbits 

and 

single -particle 

of angular 

onto spherical basis, C* 

Hartree-Fock 

energy, 

momentum 



- 

energy 

e A> 

onto z-axis, 

ld 5/2 

2s l/2 

ld 3/2 


MeV 

m A 




Ne 20 

-16. 42 

1/2 

0. 7155 

-0. 5731 

-0. 3996 

(Erf = -38. 95 MeV ) 






24 

Mg 

-18.66 

1/2 

0. 6997 

-0. 6400 

-0. 3175 

(E hf = -75.68 MeV) 

-12. 14 

3/2 

. 9616 

0 

-.2744 


> 20 24 

The C. 's for neon 20 and magnesium 24 (Ne and Mg ) have been obtained from a 

computer code (private communication, W. Bassichis) which solves the variational prob- 
lem by an iterative procedure. In table I, the results of these calculations are given for 


2 A Yukawa radial dependence was also investigated and yielded results almost iden- 
tical to the Gaussian dependence. 
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occupied proton (or neutron) orbits with positive m-projection. The orbits for the time- 
reversed states may be obtained by application of equation (3). The degenerate orbits are 
specified by the HF single -particle energies, e. (see eq. (33)). 

a 20 24 

It is clear from examination of table I that the orbits in both Ne and Mg are 
strongly deformed. These results are in agreement with results obtained by other re- 
searchers (ref. 16). 


Direct Interaction 

The calculation of microscopic form factors for inelastic scattering requires a know- 
ledge of the projectile-nucleon interaction. At the present stage of development, the 
theories of the nucleus and reactions are not sufficient to provide such a force. Instead, 
either a phenomenological interaction or the two-nucleon t-matrix must be used. From 
this point on, the projectile is assumed to be a nucleon. 

We have chosen a phenomenological spin -dependent interaction which was used by 
Glendenning and Veneroni in a recent study of inelastic scattering (ref. 4). These authors 
assumed the two-body interaction to be important only in even states, and used an inter- 
action of the form 


V(r) = -52e 


-(r/1 


95) 2 (p 


te + °- 6 P SE ) 


(60) 


where Prp E an( i Pg E are projection operators for the triplet-even and single t-even 
states, respectively. In terms of single -particle operators, 


V(r) 


-52 


0.3 1 


1 

2 





e -(r/1.85) 2 


(61) 


where Tq and refer to the isospin of the projectile and bound nucleon, respectively. 
Comparing equation (61) with equation (13), it is shown that (with r = |tq - r^|) 


= - 15 - 6 ( 1 - \ V f l)e" (r/1 ‘ 85r 


(61a) 


and 


t i( ? o» r i ) 




e -(r/1.85) 2 


(61b) 
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Because the projectile is in a pure isospin state and because the charge exchange is not 
considered an inelastic scattering event, the matrix element of • Tq will be +1 when 
projectile and target nucleon are like nucleons and -1 when they are unlike nucleons. 
Hence, the form factor can be written as the sum of the two terms, one for protons and 
one for neutrons. 


Form Factors for Neon 20 and Magnesium 24 


We have computed the form factors f^ OT for excitation of the first 2 + state of Ne 2 ® 

24 J-iocJ 

and Mg , assuming the projectile to be a proton. In figures 1 to 4, the Hartree-Fock 
form factors are compared with those obtained from the shell model and the macroscopic 

collective model. The shell-model states are obtained by filling the Id,- /0 spherical or- 

20 24 2 ^ 

bits to obtain either Ne or Mg , and then states of good J and J^. are projected 

from the resulting intrinsic state. The macroscopic form factor is proportional to the 
derivative of the real part of the optical potential which is used to fit the elastic scatter- 
ing (ref. 17). In keeping with convention, this form factor is assumed spin dependent, 
and has a derivative Woods-Saxon radial dependence: 


f • LL _ p 

where Cj depends on multipolarity of the transition, the optical model parameters, and 
the deformation /3j of the nucleus. The well parameters, Rg and ag, refer to the ra- 
dius and diffuseness of the real part of the optical potential V Q p( r o) and Tq is the sepa- 
ration between the projectile and the nuclear surface. If Uq as the strength of the opti- 
cal potential, 


Rg\ / Rg ’ 

2 + exp (r Q - — ) + exp |-r + ° 


-i 


(62) 


C J = 


^J U 0 R S 
a g (2J + 1) 1/2 


20 24 

In figures 1 and 2, the HF form factors for Ne and Mg are compared with the 

20 24 

corresponding shell-model form factors. For both Ne and Mg , the scalar parts 

(S = 0) of the HF form factors are enhanced over the shell-model results by a factor of 

about 1. 5. At small values of r(<lf) the form factors for Ne are essentially identical. 

However, as r increases, the difference becomes quite significant, and, at roughly 

r = 4f and beyond, the HF form factor is larger by approximately a factor of 2. A simi- 

24 

lar behavior is observed for the case of Mg , except that the scalar part of the HF form 
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Figure 1. - Form factors for inelastic proton scattering from neon 20. 
(The s ■ 1 curve has been magnified to make it visible. ) 
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Radius, rg, F 


Figure 2. - Form factors for inelastic proton scattering from 
magnesium 24. 


factor and the shell-model results are the same only at r < 0. 4 fermi (or fm). Another, 

> 

perhaps less important, feature is that the HF (S = 0) and the shell-model form factors 
I 20 24 

i peak at different values of r for both Ne and Mg . In both cases, the HF results 

; have their maxima nearer the accepted values of the nuclear radius (R = 1. 25 fm) 

than do the shell-model results. 

; The HF vector form factor (S = 1), which gives rise to spin-flip transitions, is 

j difficult to analyze because its shell-model counterpart vanishes identically. (This is 

23 
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due to the fact that gf^ T which appears in equation (21) is zero for the shell model 

20 24 

states. ) However, note that the HF vector form factors for Ne and Mg are about 

300 times smaller than the scalar parts. If the direct interaction which has been chosen 

(eq. (60)) is realistic, then the meaning of this result is simply that, for excitation of the 
, 20 24 

2 j states in Ne and Mg , the spin-flip part of the reaction mechanism is unimportant. 
A more reliable investigation of this question must be relegated to later studies in which 
angular distributions are compared with experiment. 

The macroscopic form factors are compared with the HF results in figures 3 and 4. 
Only the scalar part of the HF form factor is plotted, because the macroscopic model 
does not yield a vector form factor in the approximation which leads to equation (56). 

The normalizations Cj for the macroscopic form factors are based on the parameters 
which appear in the figure captions. These parameters are obtained from elastic scatter- 
ing studies (ref. 18) and the measurement of the lifetime of the 2 + states (ref. 19). 



Figure 3. - Form factors for inelastic proton scattering from neon 20. Radius of 
Wood-Saxon well, 3. 393 femtometers-, deformation parameter, 0.87; diffuseness 
of Wood-Saxon well, 0.6 femtometer; strength of real part of optical potential, 

45 MeV. 
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Figure 4. - Form factors for inelastic proton scattering from magnesium 24. 

Radius of Wood-Saxon well, 3.605 femtonete'-S; de ! ormation parameter, 

0.65; diffuseness of Wood-Saxon well, 0. 6 femtometer; strength of real 
part of optical potential, 45 MeV. 

20 24 

Figures 3 and 4 show that, although the collective form factors for Ne and Mg 
have the same shapes as the HF results, they are enhanced by factors of about 5 (for 
Ne^) and 4 (for Mg^). When the differential cross section is calculated, because it is 
related to the square of the form factor, the enhancement in magnitude of the cross sec- 
tions will be of the order of 20 times larger for the collective model. In figure 3, the 
maxima of the form factors occur at about the same radius, whereas, in figure 4, the 
macroscopic result has its maximum at a slightly larger radius. 


SUMMARY OF RESULTS 


The inelastic scattering of nucleons and composite particles from nuclei has been 
formulated in terms of a projectile -nucleon interaction and a many -body variational wave 
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function which is obtained from structure calculations. The actual nuclear states are 
generated from the variational wave function by projection - a method which has been 

quite successful in the study of low-lying levels of s-d shell nuclei. 

20 24 

The form factors for inelastic scattering of protons from Ne and Mg have been 
computed, based on the assumptions given above, and compared with the form factors 
based on the projected shell model and the macroscopic collective model. The vector 
part of the HF form factor is very small, as one might expect because the shell model 
counterpart vanishes identically. The scalar part is larger than the corresponding shell- 
model form factor. The macroscopic form factor is even larger than the HF result, 
however, indicating that cross sections which are calculated using the HF method with 
projection, will be too small by a factor of about 20. 

The model which has been used in these calculations is based on the assumption that 
the core is spherical and inert. Calculations by the present authors of B(E2) rates have 
yielded results which are also smaller than experiment. This is in agreement with the 
form factor study presented for both Gaussian and Yukawa interactions. Thus, the HF 
theory appears to be inadequate to provide the correct magnitude for cross sections, 
etc. , if the inert core assumptions is retained. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 9, 1968, 

129-02-07-07-22. 
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